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Abstract 

A new time projection chamber (TPC) was developed for neutron lifetime 
measurement using a pulsed cold neutron spallation source at the Japan Pro¬ 
ton Accelerator Research Complex (J-PARC). Managing considerable back¬ 
ground events from natural sources and the beam radioactivity is a chal¬ 
lenging aspect of this measurement. To overcome this problem, the devel- 
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oped TPC has unprecedented features such as the use of polyether-ether- 
ketone plates in the support structure and internal surfaces covered with 
6 Li-enriched tiles to absorb outlier neutrons. In this paper, the design and 
performance of the new TPC are reported in detail. 
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1. Introduction 


Neutron lifetime is an important observable parameter used to determine 
V u d in the Cabibbo-Kobayashi-Maskawa (CKM) matrix, together with the [3- 
asymmetry parameter in neutron decay. Furthermore, it is used as a probe to 
test the Big Bang theory through primordial nucleosynthesis. After the first 
observation of neutron decay in 1948, various neutron lifetime measurements 
using nuclear reactors have been conducted. Currently, two main measure¬ 
ment methods exist with a 8.4 ±2.2 s discrepancy between their results: one 
requires the counting of surviving ultra-cold neutrons after storing (giving 
879.6 ± 0.8 s [l|, 2, 3, 4, 5]) and the other requires the counting of trapped 
protons from neutron decay (giving 888.0 ±2.1 s {s, 7]). Therefore, accurate 
measurements using different methods are necessary. 

The proposed experiment in this paper firstly employs an accelerator 
as a neutron source [8|, which is conceptually based on a measurement by 
Kossakowski et al. jsj] using a reactor and a time projection chamber (TPC). 
The result of this measurement was published in 1987 as 878 ± 27 (stat.) ± 
14 (sys.) s. In this method, the TPC is filled with 4 He, CO 2 , and a few 
ppm of 1 * 3 He gas, and detects electrons emitted from the neutron decay while 
simultaneously measuring the neutron flux by counting 3 He(n, p) 3 H reactions. 
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To keep the number density of the 3 He, the TPC is housed inside a vacuum 
vessel, which is filled with gas after vacuuming of the vessel and sealed during 
operation. 

The incident neutron beams are shaped into short bunches with lengths 
of approximately half the TPC. The fiducial time during which the neutron 
bunch is entirely inside the TPC can be defined, which enables us to reduce 
the uncertainty related to comparison between the number of the neutron 
decays and the number of the 3 He(n, p) 3 H reactions. 

2. Experimental overview 

In this section, our experiment is described in comparison with the exper¬ 
iment conducted by Kossakowski et al. The specifications of the beam and 
the TPC are summarized in Table [0 In particular, the TPC developed in 
this paper had two features which can be found in Table [2] and [3] in Section 

El 


• made of PEEK without radioactive contamination, 

• lined by 6 Li tiles for the capture of scattered neutron. 

Our TPC achieved the drift velocity of 1.0cm//rs under 100 kPa. The multi¬ 
plication gain was 4 x 10 4 with a resolution defined by FWHM of 22.9% for 
5.9 keV, which resulted in 0.2 keV energy threshold per wire. The detailed 
performance is described in Section 0 

The coordinate system is defined as follows: the z-axis is parallel to the 
direction of the neutron beam in the TPC, the y-axis is vertically aligned 
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from the bottom to the top, and the x-axis is defined using the right-handed 
Cartesian coordinate system. 


2.1. The Kossakowski et al. neutron lifetime experiment 

Kossakowski et al. employed a continuous cold neutron beam from a 
reactor at the Institut Laue-Langevin (ILL), with a TPC of volume 190 
(x) x 190 (y) x 700 (z) mm 3 10]. The neutron beam was shaped into 
bunches with lengths of 23-25 cm by a chopper drum rotating at 110 Hz, 
and monochromatized at a wavelength of 4.73 A by Bragg reflection on a 
graphite crystal llj. The size and the divergence of the neutron beam were 
15 x 25 mm 2 and ±8.7 mrad, respectively. The neutron flux inside the TPC 
was 2.2 x 10 5 /s, which corresponded to 0.10 neutron decays/s for the fiducial 
time of 400 /is. A duty factor was calculated as 0.044. 

The TPC consisted of a drift cage with a multi-wire proportional chamber 
(MWPC) inside the vacuum vessel. The total gas pressure was fixed at 95 
kPa, and a mixture of 93% 4 He, 7% C0 2 and 0.7 ppm 3 He was adopted. The 
MWPC had sense/held wires in the z-direction sandwiched by two layers of 
cathode wires in the x-direction. The wire cell in the MWPC was 10 x 10 
mm 2 . The data acquisition was triggered by any hit on the sense wires. 

Electrons from the neutron decays have a continuous kinetic energy spec¬ 
trum up to 782 keV and deposit a part of this energy in the TPC; The energy 
loss for the electrons with kinetic energy of 0(100) keV is less than 1 keV/cm. 
On the other hand, the 3 He(n, p) 3 H reactions release monochromatic Q-value 
energy of 762 keV, and its decay products are both stopped inside the TPC. 
Due to the saturation of the multiplication at the sense wires, as described 
in Appendix A the energy spectrum of the 3 He(n, p) 3 H reaction becomes 
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broad and overlaps that of the neutron decay. Kossakowski et al. used the 
maximum pulse amplitude among the sense wires as a discriminant vari¬ 
able and set 120 keV for the energy threshold. An uncertainty of 0.6% was 
assigned for the separation. 

The neutrons are absorbed by the carbon in the gas at a rate comparable 
to the neutron decay. The 12 C(n, 7 ) 13 C reaction generates a point-like energy 
deposit of 1.0 keV along the beam axis due to the 13 C recoil by prompt 7 - 
ray of 5.0 MeV. To prevent these events, more than two sense wire hits 
were required. Thus, e car bon represents the fraction of the loss due to the 
removal of the 12 C(n, y) 13 C reactions, which was calculated as 11 ± 1% in the 
Kossakowski et al. experiment. 

In this approach, the triggered event data were divided into two subsets: 
lV 3 He for the 3 He(n, p) 3 H reactions and Np for the neutron decays. The de¬ 
tection efficiency for the neutron decay, ep, was more than 99.9%, so that the 
uncertainty on ep was less than 0.1%. The neutron lifetime, r n , is expressed 
as 


^ £/3 X (1 ^ backscatter ) X (1 6carbon ) 

Np B en y -Bfteam P 3 He^ 3 He^n 

where /? 3 He , u 3 He , and v n are the number density of 3 He, the cross section 
of the 3 He(n, p) 3 H reaction, and the neutron velocity, respectively. During 
filling the gas, the value of p 3 He was determined by a pressure gauge with 
an uncertainty of 0.5%, while cr 3 He u n was treated as a constant using the 
1/v dependence of < 73 He for the cold neutron. At v n = 2,200 m/s, cr 3 He was 


obtained as 5, 333 ± 7 barn 


12]. The remaining terms are described below. 


Neutron decay events were identified only by detecting electrons by the 
TPC, thus, there are many types of natural environmental background events 
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labeled B env : cosmic rays, 7 -rays from natural radioisotopes, and radioactiv¬ 
ity due to the detector materials. Their TPC with cosmic-ray veto counters 
and lead shields had B env of 80 cps, since no special low-activity materials 
were used. 

I 11 Eq.flT]) above, Ng includes background events labeled B beam due to 
prompt 7 -rays induced by the neutron beam, which occured at 20 cps. The 
7 -rays from the upstream side of the TPC were prevented by the lead shields. 
On the other hand, once the neutrons passing through the TPC were scat¬ 
tered by the gas and absorbed in the TPC itself or the vacuum vessel, they 
also became a 7 -ray source. The TPC was surrounded by a 6 LiF coated 
plexiglas box, since 6 Li has a large neutron absorption cross section but does 
not emit prompt 7 -rays. However, neutron absorption by the TPC structure 
itself was not avoided in this setup. 

I 11 order to reduce this background events, one of the terminal points of 
the reconstructed tracks was required to be located on the neutron beam 
axis, since the electrons emitted by the neutron decays originate there. This 
requirement obstructed some of the electrons originating from the neutron 
decay which were scattered on the inner surface of the TPC. The term, 
ebackscatter, represents the fraction of the loss, which was estimated by Kos- 
sakowski et al. as 1.5 ± 0.5%. 

Also, the large B env limited the statistical uncertainty to 3.1%. Kos- 
sakowski et al. then estimated the value of B beam as 0.9% of the neutron 
decay rate and added the same amount of systematic uncertainty to the sub¬ 
traction of B beam from the result. All the uncertainties for the measurement 
are summarised in Table [4j 
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2.2. Neutron beam for J-PARC experiment 

The Materials and Life Science Facility (MLF) at the Japan Proton Ac¬ 
celerator Research Complex (J-PARC) has one of the most intense pnlsed 
neutron sources currently available. 3 GeV proton beams with a repetition 
rate of 25 Hz are injected into a liquid mercury target, which is supported 
from the downstream of the proton beam line. The spalled neutrons emitted 
from the target are cooled by liquid hydrogen moderators located above and 
below the target. Our experimental setup, located at the polarized-beam 
branch of the neutron optics and physics (NOP) beamline BL05 port 13], is 
shown in Figure [0 

The polarized pulsed neutron beam has a broad time-of-flight (TOF) 
structure reflecting Boltzmann distribution. The polarized beam is trans¬ 
ported into the TPC via a neutron optical system of magnetic supermirrors 
and spin-flipper coils called a spin flip chopper (SFC) [14]. Resonance spin 
flipper coils can rotate the spin orientation of the neutrons, and the flipped 
neutrons then pass through the supermirrors and are dumped before the 
TPC. By turning on the coils during the un-desired time domain, the beam 
is shaped into arbitrary-sized bunches. The neutron flux is monitored using 
a thin 3 He beam monitor 15] at the end of the SFC. The beam size is finally 


defined as 20 x 20 mm 2 at the neutron switching shutter positioned between 
the SFC and the TPC. 

Assuming a proton-beam power of 300 kW at J-PARC, the number of 
neutrons at the exit of the polarized-beam branch was calculated to be 2.9 x 
10 7 /s with an aperture of 30 x 30 mm 2 . The divergence of the beam was 
±4.2 mrad. If the entire TOF region passed the SFC, 1.2 x 10 6 /s neutrons 
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Figure 1: Experimental setup at BL05 in MLF (top view). (A) BL05 concrete shield, 
(B) lead shield, (C) iron shield, (D) 6 Li beam collimator, (X) polarized-beam branch, 
(Y) unpolarized-beam branch, (Z) low-divergence-beam branch, (a) short-wavelength pass 
filter, (b) guide coil, (c) resonance spin flipper coils, (d) magnetic supermirrors, (e) neutron 
beam monitor, (1) 50-/mi-thick Zr window, (2) neutron switching shutter, (3) cosmic-ray 
veto counter, (4) lead shield, (5) vacuum vessel, (6) TPC, (7) 6 Li beam catcher, and (8) 
turbo molecular pump. 






















































Figure 2: Schematic view of our TPC. The TPC support structure is composed of PEEK 
plates. The 6 Li tiles are set inside the drift cage and above the MWPC. A neutron bunch 
have a length of approximately half that of the TPC. 

were guided inside the TPC. When the SFC reflected five 40-cm-long bunches 
per pulsed beam, 1.7 x 10 5 /s neutrons were obtained. In this configuration, 
the duty factor was calculated as 0.059. 

2.3. TPC for the J-PARC experiment 

The vacuum vessel was made from a stainless steel housing with internal 
size of 630 (x) x 680 (y) x 1180 (z) mm 3 and three aluminum lids. The 
TPC comprised a MWPC on a drift cage as shown in Figure [2j The volume 
of the drift cage was 290 (x) x 295 (y) x 960 (z) mm 3 , and 0.092 neutron 
decays/s were expected in the fiducial time during which the neutron bunch 
was entirely inside a region of ±400 mm at the center of the TPC. The wire 
structure of the MWPC was the same as that of the TPC developed by 
Kossakowski et al, and the wire cell size in the MWPC was 12 x 12 mm 2 . 
The expected neutron decay rate of 0.092 cps in the fiducial time was 
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equivalent to the continuous events of 1.6 cps. For further precision regarding 
the statistical uncertainty, reduction of B env was crucial, especially in relation 
to the background events due to radioisotopes contained in the TPC system. 
We investigated several candidates for the TPC material and decided to use 
PEEK for the main structure as described in Section 13.11 

In order to reduce the systematic uncertainty concerning the subtraction 
of B beam , we developed a new tile containing 6 Li as shown in Section 13721 and 
covered the entire inner surface of the TPC with these tiles. The pressure 
of the gas inside the TPC can be varied from 50 to 100 kPa in order to 
evaluate B beam using a data-driven method. If we use different pressures for 
the measurement, p\ and p 2 , the neutron lifetime can be expressed by Eq. tfTj) ; 


= 


^SHeOl 


X 


e(Pi) 


B env (p | ) Bffgamipi) P 3 He^" 3 HeKi 


^3He(P2) 


X 


e(P2 


( 2 ) 


(P‘2 ) B env (p 2 ) B beam iP'2 ) P 3 Hc*^" 3 He^n 

where e denotes ep x (1 - e backscatt er) x (1 -e carbon ) for simplicity. Since B beam 
is proportional to the pressure, we have 

BbeamiPl) _ Pi 

B()enni iP'l ) P : 2 

From Eq.(j2]) and Eq.(j3J), the neutron lifetime can be calculated without B beam 
as 

Pie(p 2 )N 3Uc (p 2 ) - P 2 e(pi)N 3lie (p 1 ) 1 


r„ = 


x 


( 4 ) 


Pl{Np(p 2 ) - B env (p 2 )) - P2(Np(pi) - B env (p 1 )) pSHeCTSHeUn’ 
which would reduce the uncertainty related to B beam . In this study, the gas 
mixture was fixed at He/C02 = 85/15. The 3 He abundance was controlled 
as about 1 ppm. 
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Reduction of the systematic uncertainty regarding the correction for B carbon 
is possible by introducing an energy threshold to distinguish the 12 C(n, 7 ) 13 C 
reaction from the neutron decay, instead of just requiring more than two sense 
wire hits. 99.9% of electrons from the neutron decay have kinematic energy 
of more than 4.0 keV. If an energy threshold of 4.0 keV would be applied, 
the B car b on correction would be reduced to 0.1%, so that the uncertainty on 
this correction would also be reduced to less than 0.1%. The detailed design 
and performance of the TPC are discussed in Sections 0] and [HI respectively. 


3. TPC material selection 

The PEEK and the 6 Li tiles are the unique materials used for the TPC de¬ 
sign presented in this paper. In this section, their specifications are described 
in detail. 


3.1. PEEK as the TPC support structure 

Activity concentrations of the radioisotopes in the materials were deter¬ 
mined by y-ray spectrometry using a germanium detector in KEK, which 


was previously calibrated by Canberra 


16] and shielded from environmental 


radiations by a 0.5-cm-thick acrylic plate, 0.5-cm-thick oxygen-free copper, 
and 10-cm-thick lead. The 7 -ray count rate from the sample material was 
obtained by subtracting the count rate measured without the material. The 
detection efficiency of each 7 -ray was determined by using Canberra’s Lab- 


17]. The nuclear data obtained from the Table of Isotopes 


SOCS software 

Q were used for data analysis. 

Table [2] shows the properties of Polyether-ether-ketone (PEEK) and Poly¬ 
phenylene-sulfide (PPS) produced by Yasojima Proceed Co. Ltd with product 
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Kossakowski et al. 

This work 

Facility 

ILL 

J-PARC (300 kW) 

Repetition rate [Hz] 

110 

25 

Pulse per repetition 

1 

5 

Beam size [mm 2 ] 

15 x 25 

20 x 20 

Beam Pulse length [cm] 

23-25 

40 

Beam divergence [mrad] 

±8.7 

±4.2 

Velocity [m/s] 

837 

500-1200 

Duty factor for fiducial time ( —F ) 

0.044 

0.059 

Neutron flux inside the TPC [1/s] 

2.2 x 10 5 

1.7 x 10 5 

Neutron decay rate [l/(F-s)] 

0.10 

0.092 

Drift cage size [mm 3 ] 

190 x 190 x 700 

290 x 295 x 960 

MWPC cell size [mm 2 ] 

10 x 10 

12 x 12 

TPC Gas pressure [kPa] 

95 

50 ~ 100 

Gas mixture ratio [He : C0 2 ] 

93 : 7 

85 : 15 

3 He abundance [ppm] 

0.7 

~ 1 


Table 1: Comparison of neutron beams and TPCs between the experiment by Kos- 
sakowski et al. and the present experiment, assuming 300-kW beam power at J-PARC. 
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names PEEK450G and FORTRON, which had the Erst and the second least 
activity concentration emitting 7 -rays among the candidates for the TPC 
support structure. The PEEK possessed enough mechanical properties in 
terms of the wire tension, soldering and vacuuming. Table [2] also shows the 
activity concentration of the 6 Li tile described in Section 13.21 stainless steel 
(SUS304) and aluminum (A5052) for the vacuum vessel provided by VIC 
International, Inc. 

The dominant radionuclides inducing the background events in the TPC 
were 210 Pb and its daughters, 210 Bi and 210 Po, according to the 7 -ray spec¬ 
trometry. The radiations (a-rays, /3-rays, 7 -rays, and bremsstrahlung) emit¬ 
ted following the decay of 210 Pb and its daughters in the PEEK and the 
SUS304 were simulated by the Geant4 Q with the measured activity con¬ 
centrations and the detailed geometry of the detector system as described 
in Section [4j The activities of the daughter radionuclides were assumed to 
be in equilibrium with 210 Pb. The trigger rates of the TPC induced by the 
radiations from 210 Pb and its daughters were anticipated to be 2.2 ± 1.5 cps 
from the TPC support structure and 0.2 ± 0.1 cps from the vacuum vessel. 
The uncertainty was derived from the measurements of the activity concen¬ 
tration by the germanium detector. I 11 total, the expected background rate 
was 2.4 ± 1.5 cps. Comparison with the observation can be shown in Section 
I5R1 

3.2. 6 Li tiles inside the TPC 

To produce the 6 Li tile, 6 Li-enriched lithium fluoride (LiF) powder was 
baked with poly-tetra-fluoro-ethylcne (PTFE). The 6 Li had a 95 % concen¬ 
tration in the lithium. The mixture ratio of the LiF was 30% by weight 
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Total 7 -radio- 

nuclides 

[Bq/cm 3 ] 

210 Pb 

[Bq/cm 3 ] 

Elastic 

modulus 

[GPa] 

Melting 

point 

PC] 

Water 

absorption 

m 

PEEK 

0.017 ±0.012 

0.022 ± 0.011 

3.6 

334 

0.14 

PPS 

0.057 ±0.015 

0.053 ±0.015 

3.9 

278 

0.04 

6 Li tile 

- 0.011 ± 0.022 

- 0.010 ± 0.022 

— 

— 

— 

SUS304 

1.5 ±1.0 

1.5 ±1.0 

— 

— 

— 

A5052 

-0.019 ±0.080 

-0.028 ±0.080 

— 

— 

— 


Table 2: Total activity concentration of the radionuclides emitting 7 -rays and activity 
concentration of 210 Pb in the materials for the TPC and the vacuum vessel, determined 
by 7 -ray spectrometry. Mechanical properties of PEEK and PPS are taken from the 
specification sheet of PEEK450 and FORTRON produced by Yasojima Proceed Co. Ltd. 


due to the stiffness property. The thickness of the 6 Li tile was 5 mm, and 
its molding size before manufacturing was 300 x 300 mm 2 . The absorption 
length of the 6 Li tile was 0.5 mm for the thermal neutrons, so that the neu¬ 
trons scattered by the gas can be completely absorbed. As Table [3] shows, 
the ratio of the neutron absorption with the prompt 7 -rays in the 6 Li tile 
was reduced to 8.3 x 10 -5 calculated by using cross sections in ref. 20]. If 
the prompt 7 -rays are emitted, the average number of the prompt 7 -rays was 
calculated as 2.0, which was based on the NNDC on-line data service from 


the ENSDF database 2lf]. In the developed TPC, 100 -/irri-thick sheets of 
PTFE were attached to the 6 Li tile in order to prevent a and 3 H particles 


from returning to the TPC active volume. 
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Absorption cross 

section with 7 -ray 

Mole 

fraction 

Branching 

ratio 

average number 

of 7 -ray 

6 Li 

39 mbarn 

0.17 

4.1 x 10 “ 5 

1.4 

7 Li 

45 mbarn 

0.01 

2.5 x 10 “ 6 

1.0 

C 

3.5 mbarn 

0.21 

4.5 x 10 " 6 

1.2 

F 

39 mbarn 

0.61 

3.5 x 10 " 5 

2.8 

6 Li tile 

32 mbarn 

1.0 

8.3 x 10 " 5 

2.0 


Table 3: Properties of the prompt 7 -rays induced by neutron capture in the 6 Li tile. 

4. The detector system 

In this section, the structure of the drift cage and the MWPC in our TPC 
are explained in detail. Figure [3] shows a cross section of the TPC. 

4-1. Drift cage 

The support structure of the drift cage was composed of PEEK plates 
which were welded into frames for the top and base sides, horizontal planks 
for the longer side, and perpendicular panels with a rectangular neutron- 
beam-passing hole for the shorter side. At the corners of the cage, 5-mm- 
pitch drift wires (050-/im of BeCu alloy) were soldered on a polyimide sheet 
with an etched electric circuit pattern. 2-Mfi chip resistors interconnected 
the drift wires. 

The 6 Li tiles were fixed inside the PEEK structure. The base and longer 
side surfaces were covered with three and half 6 Li tiles, and the shorter side 
surfaces were covered with one 6 Li tile. The edges of the 6 Li tiles were milled 
so that the adjacent 6 Li tiles overlaped by 5 mm in order to close any gaps. 
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Figure 3: TPC in vacuum vessel viewed from downstream side (x-y plane). 

A 12-/jm-thick aluminized PET film was added to the surface of the bottom 
6 Li tiles, for the application of a negative voltage for drifting of the ionized 
electrons. 

The presence of dielectric materials such as the 6 Li tiles and the PEEK 
structure inside the drift wires could distort the drift field, and this effect 
must be considered in the design of the drift cage. The relative permit¬ 
tivities of the 6 Li tile and the PEEK structure were measured as 3.0 and 
3.2, respectively. We then simulated the drift field using a three-dimensional 
static electric field analysis to make the equipotential surface parallel to xz- 
plane. The calculated drift voltage was uniform within 1.7% in the cage. 
Note that wire structures of MWPC and charging up of materials were not 
taken into account for this simulation. 

A 40 x 40-mm 2 and a 60 x 60-mm 2 apertures were placed at the entrance 
and exit of the cage, respectively, to allow the neutron beam to pass, taking 


16 




































the divergence of the beam into consideration. The beam ducts were con¬ 
nected to them, and the 6 Li tiles covered the inner surfaces of the ducts as 
well as the TPC. As a beam catcher, a PTFE box filled with the 6 Li-enriched 
LiF powder was located at the end of the beam duct. 

J h 2. The MWPC 

The three-layer MWPC wire structure with 6-mrn gaps was made on two 
6 -mm-thick PEEK frames. The outer and inner sizes of the PEEK frames 
were 412 (x) x 1,070 (2) mm 2 and 312 ( x ) x 970 (z) mm 2 , respectively. A 
single frame was used to stretch both the sense/hcld wires in the ^-direction 
and the top cathode wires in the ^-direction; the other frame had the bottom 
cathode wires on one side. An additional PEEK frame was mounted above 
the MWPC so that the top 6 Li tiles could be fixed upon it. The top 6 Li 
tiles were covered with aluminized PET him, for the application of a positive 
voltage, such that the ionized electrons above the MWPC inside the TPC 
were collected into the him. 

In the MWPC anode layer, 28 sense wires (020-/mi of gold-plated tung¬ 
sten) and 27 held wires (050-fim of BeCu alloy) were alternatively placed 
at every 6-nnn pitch in the ^-direction. Each cathode layer consisted of 160 
6 -mm-pitch cathode wires (050-fim of BeCu alloy) in the x-direction. The 
wire positions were determined by 100-/mr-width gutters. 

The signals from 24 sense wires and 24 held wires in the central region 
were individually read out. Four sequential cathode wires were linked to 
a single channel, which corresponded to 40 channels per cathode layer. In 
total, the TPC had 128 channels. 

Charges induced on the MWPC wires were converted to the voltage signal 
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by charge-sensitive amplifiers. The amplifiers were mounted on the top of 
the TPC. Two conversion factors of 1.0 V/pC (high gain) and 0.1 V/pC (low 
gain) were implemented. The high-gain amplifiers were for the sense and 
bottom cathode wires, and the low-gain amplifiers were for the field and top 
cathode wires. It was confirmed that the 3 He(n, p) 3 H events can be recorded 
without exceeding the dynamic range, because of the space-charge effects 
occurring during the multiplication. 


4-3. Data acquisition system 

A common pipelined platform for electronics readout (COPPER) devel¬ 


oped at KEK 


22] was employed for our experiment. We used COPPER lite, 


which was a successor to the COPPER device. A COPPER board can mount 
up to four digitization daughter cards, called front-end instrumentation en¬ 
tities for sub-detector specific electronics (FINESSE). 

The signals from the 24 TPC sense wires were fed into discriminators with 
20-mV threshold for the pulse height, and a logical sum of the 24 outputs 
from the discriminators triggered the data acquisition. A TDC with a 1.25- 
MHz clock and 16-bit dynamic range was employed to measure the time 
period from the injection of the latest proton pulse on the mercury target 
and the generation of the trigger signal. The waveforms of the 128 channels 
from the TPC were digitized using a flash-ADC-type FINESSE card with 
a 20-MHz clock, 12-bit dynamic range, and 512-words/channel FIFO. 


4-4■ Calibration system 

To monitor the TPC conditions, a 55 Fe source with a collimator was 
equipped on a vertical rotation stage as shown in Figure El Two slits at 
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Figure 4: Box diagram of the trigger and data acquisition. 

the side of the drift cage can be selected for the injection of 5.9 keV X-rays, 
with distances from the MWPC of 75 and 225 mm. The X-rays passing 
through the collimator spreaded in the xz-plane in a trapezoidal shape, with 
a ^-divergence of 10 mrad and covering 60% of the MWPC area. When the 
source was not used, the slits were closed by small b Li tiles that were also 
located on the rotation stage. 

4-5. Shields 

The vacuum vessel was housed in a lead shield which covered all six sides 
except the beam path. The lead was of 10-cm thickness on the upstream 
side, and 5-cm-thickness on the other five sides. The hermeticity of the lead 
shield was 99 % because of the holes for the beam ducts. The lead shield 
was further reinforced by a 20-cm-thick iron shield on the upstream side and 
on the side of the TPC that met another beam line, “BL04”. According 
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to the Nal detector measurement, 98 % of the environmental 7-rays were 
suppressed. 

4-6. Cosmic-ray veto counters 

The lead shield surfaces, apart from the bottom, were covered with cosmic- 
ray veto counters. Both the lead shield and the cosmic-ray veto counters had 
square holes around the beam duct. A panel of cosmic-ray veto counters 
consisted of a two-layer array of extruded plastic scintillator boards, with 
the layers being optically isolated. The scintillator had trenches with 20-mm 
pitch, and a wavelength shifting optical fiber was inserted in each trench. The 
fibers of the layers were bundled and read out by a photomultiplier tube. 

The performance of the cosmic-ray veto counters was evaluated using a 
cubic scintillator with a side length of 10 cm, which was also placed on the 
center of the bottom lead shield. The cosmic rays deposited approximately 20 
MeV in the cubic scintillator, while other environmental background sources, 
including the neutron-induced background, deposited at most 8 MeV. Thus, 
pure cosmic rays could be extracted. The measured efficiency of the cosmic- 
ray veto counters was found to be 96%. 

4- 7. Apparatus for determination o/ 3 He partial pressure 

To determine the 3 He partial pressure in the TPC, we used a pressure 
gauge produced by Mensor Corp. (Mensor Digital Pressure Gauge Series 
2500), which was equipped with a piezoresistive transducer. This gauge 
achieved 0.01% accuracy for atmospheric pressure, as confirmed by NIST- 
traceable standards. Current precision for the determination of the 3 He 
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abundance of 1 ppm is 0.3%. We have also developed a sector-type single- 


23] in order to directly 


focusing mass spectrometer in the University of Tokyo 
measure the 3 He and 4 He ratio. The mass spectrometer employed a double 


collector system for simultaneous detection of 3 He and 4 He. 


5. TPC performance 

In this section, we discuss in detail the detection system performance in 
terms of various TPC characteristics. 

5.1. Detector configuration 

For measurement of the drift velocity, the multiplication gain, and the 
detection efficiency, the top 6 Li tiles on the TPC were removed. At that 
time, the MWPC was covered with ethylene vinyl alcohol (EVOH) film to 
obstruct ionized electrons from above the TPC, so that the particle response 
of the TPC would be unaffected. For the evaluation of the long-term stability 
and background event rate, the top 6 Li tiles on the MWPC were mounted in 
the normal setup. 

According to Eq. (ID . the measurements of the neutron lifetime under 
several gas pressure could reduce the systematic uncertainty on the subtrac¬ 
tion of B beam . We investigated the detection efficiency under 50 kPa and 
100 kPa. For the drift velocity and the multiplication gain, 75 kPa was also 
evaluated. 

5.2. Drift velocity 

Since the total count rate of the TPC was about 1 kcps which was dom¬ 
inated by the cosmic rays, the memory time of the TPC should be less than 
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a few tens /zs in order to keep the dead time of the data acquisition below 
0 (1)%, and the drift velocity should be faster than about 1 cm f[is with the 
maximum drift distance of 295 mm of the TPC. 

The drift velocities of the ionized electrons in the TPC, Vd, were measured 
using the cosmic rays which traversed the TPC in the ^-direction. The timing 
difference, At, between the earliest and the latest hits on the sense wires was 
calculated. At enc i was defined as the end point in the distribution of At. 
Due to the pitch of the anode wire, the apparent drift distance, Vd x At end , 
was gradually close to 295 mm, if the angle between the cosmic ray and the 
drift direction would be rectangular. Thus we analyzed the cosmic rays with 
different angles and took the limit for the horizontal cosmic ray in order 
to obtain We implemented this method in the simulation, and checked 
an output of the calculation for a known velocity of the cosmic ray as an 
input. Systematic uncertainty of Vd was set as 4%, which was determined by 
taking the maximum deviation between the inputs and the outputs of the 
simulation. 

Figure [5] shows Vd as a function of the reduced electric held for the gas 
configuration, which is compared with the Magboltz simulation 24]. Tem¬ 
perature in the vacuum chamber varied from 298 K to 309 K during the 
measurements, thus we draw the results as a band calculated with the mini¬ 
mum and maximum temperatures. Under 50 kPa and 100 kPa gas pressure 
for the other measurements in this paper, the high voltage for the drift cage 
was set as —6, 000 V and —9, 000 V to satisfy the requirement, which cor¬ 
responded to Vd of 1.4 cm / fxs and 1.0 cm//is, respectively. The measured 
velocities were slightly higher than Magboltz simulation, but could not affect 
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Reduced Drift Field E/p (V/cm/kPa) 


Figure 5: Drift velocity dependence on reduced drift field and gas pressures, where E 
and p are the electric field and total pressure, respectively. The gas mixture was fixed at 
He/C02 = 85/15. The band shows the calculated values with the Magboltz simulation 
including an uncertainty of temperature 241. 


the measurement of the neutron lifetime, since we would use the measured 
velocities themselves. 


5.3. Multiplication gain 

The multiplication gain was measured using the energy deposit of 5.9 keV 
X-rays from the 55 Fe source. The 55 Fe source was placed directly above the 
MWPC, such that the X-rays were directly absorbed in the MWPC wire-cell 
volume without drift. A W -value of 42 eV was used as an average energy 
to produce an ionized electron in the gas. Figure [6] shows that our TPC 
achieved the multiplication gain of 4 x 10 4 , which allows us to apply the 
threshold of 0.2 keV by the signal-sensing threshold at 20 mV. The high 
voltages for the sense wires under 50 and 100 kPa were set to be 1,440 and 
1,750 V, respectively. 

Pulse height spectra for 5.9 keV X-ray on a central wire can be shown in 
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Figure [7] The background events were subtracted by rotating the 55 Fe source 
stage. By fitting 2 a regions of the peaks, the energy resolutions for 5.9 keV 
defined by the FWHM were obtained as 26.1% with 50 kPa at 1,400 V and 
22.9% with 100 kPa at 1,720 V. These resolutions would realize sufficient 
separation between the neutron decay and the 1.0-keV 12 C(n, y) 13 C reaction. 
The reduction of the uncertainty on e car bon is discussed in Section [575] together 
with the long-term stability. 


30 r 
20 


Discharge limit 


50 kPa 75 kPa 100 kPa 


He/CO, = 85/15 
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Figure 6: Multiplication gain dependence of the MWPC on voltage for sense wires and 
total pressure values. The gas mixture was fixed at He/C02 = 85/15. Discharge limit was 
empirically obtained. 


5.f. Detection efficiency for minimum ionizing particles 

The detection efficiency for the minimum ionizing particles with minimum 
path length on the sense wires was measured using the cosmic rays. At the 
normal position of the TPC, the required direction was achieved only by 
horizontal cosmic rays, but the numbers were not sufficiently high. In order 
to use perpendicular cosmic rays instead, we rotated the TPC by 90° around 
the beam axis in the vacuum vessel, which has an almost square cross section 
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Figure 7: Pulse height spectra of 55 Fe X-ray on a central anode wire taken with 1,400 V 
at 50 kPa and with 1,720 V at 100 kPa. The resolutions in FHWM were 26.1% and 22.9% 
for 50 kPa and 100 kPa, respectively. 

(see Figure[3]). After the rotation, one side faced np and the other side turned 
to the bottom, while the MWPC plane was perpendicular to a horizontal 
plane. Two plastic scintillators for triggering were set both over and under 
the vacuum vessel. Cosmic rays with a zenith angle of less than ±5° were 
then selected. The detection efficiency of the sense wire was determined by 
taking the hit ratio when both adjacent sense wires showed hits. 

The expected average energy deposit on one sense wire was found to be 
0.5 and 1.0 keV for 50 and 100 kPa, respectively, while the hit discrimination 
threshold was 0.2 keV for both cases. The results are presented in Figure 
[HI The average efficiencies in the center of the TPC were 78 % for 50 kPa 
and 97 % for 100 kPa. Using the Geant4 and implementing these efficiencies, 
the detection efficiency for the neutron decay, e^, under both total pressure 
values was estimated to be over 99.9%, so that the uncertainty on ep should 
be less than 0.1%. 
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Figure 8: Detection efficiency for the minimum ionisation particles with minimum path 
length on sense wires. The gas mixture was fixed at He/CC >2 = 85/15. 

5.5. Long-term, performance stability 

Since the gas is sealed during operation of the TPC, the attachment of 
the ionized electrons in the drift cage is caused by contaminations due to 
outgassing. The 55 Fe source on the rotation stage provides two drift lengths, 
which yields an attachment coefficient, C [m -1 ]. A large attachment coeffi¬ 
cient means degradation of the energy resolution. 

The outgassing rate was related to the pressure achieved during the evac¬ 
uation before the gas filling, P vac , and 3.6 x 10~ 4 Pa was reached after an 
evacuation of approximately 2 weeks. In this condition, the outgassing rate 
was 4.4 Pa/day. Figure [9] shows the long-term evolution of the attachment 
coefficient for the two different P vac . The TPC achieved C < 0.25 m -1 for 
5 days of continuous operation, which resulted in an energy threshold of 
1.4 keV to reject the 1.0-keV 12 C(n, 7 ) 13 C reaction. With regard to e car bcm, a 
99.9% efficiency for electrons from the neutron decay can be obtained with an 
energy threshold of 4.0 keV, so that the uncertainty of the e car b on correction 
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should be less than 0.1%. 



Figure 9: Time variation of transportation efficiency of the TPC with different achieved 
pressures. The gas mixture was fixed at He/C02 = 85/15. 

The dielectric materials such as the 6 Li tile were exposed inside the TPC, 
which may result in distortion of the electric held by charging. During the 
5-day operation, the deviation of the drift velocity and the multiplication 
gain were found to be 1% and 10%, respectively, which would not affect 
the neutron lifetime measurement, since the variation can be monitored and 
calibrated correctly. 

5.6. Background event rate 

In the MLF, the value of B env without any rejection of the environmental 
y-rays and cosmic rays was 123.7 cps. To evaluate the background events, 
the same voltage as in the detection efficiency measurement was applied. 
B env was reduced to 58.4 cps by the lead shields. Finally the cosmic ray 
counters reduced B env to 7.7 cps. The count rates of the TPC due to the 
environmental y-rays and cosmic rays would be 1.3 and 2.0 cps, which were 
calculated from the performance of the lead shield and the cosmic-ray veto 
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counters described in Section [U The remaining background rate of 4.4 cps 
would be due to the radioisotopes in the TPC and the vacuum vessel, which 
was almost consistent with 2.4 ± 1.5 cps of the expectation shown in Section 
EU With the improved B env , a statistical uncertainty of 1.0% can be achieved 
in 34 hours, which is the same data acquisition time as that of Kossakowski 
et al. 

Further, the 6 Li tiles inside the TPC can reduce Bb eam to 2.0 cps. Al¬ 
though the current systematic uncertainty on the subtraction of Bb ea m has 
not yet been measured, it should be less than 0.9% because of the significant 
reduction of Bb eam . All the uncertainties on the neutron lifetime measure¬ 
ment, including the improvements, are summarized in Table [4j 

6. Conclusion 

A new measurement approach is being used at J-PARC in Japan, which 
has several unique advantages compared to the current measurement meth¬ 
ods. The TPC with PEEK plates and 6 Li tiles was developed to enhance 
neutron decay detection during long term operation using a low background 
environment and improved signal efficiency. Our experiment with this TPC 
will be able to reach 1.6% of the precision on the neutron lifetime in 34 
hours. The expected systematic uncertainty is dominated by the subtraction 
of the background due to the neutron beam interacting with the gas inside 
the TPC, which can be improved by the comparison of the measurements 
among the different pressures of the gas. We have proved the TPC to keep 
the performance under the pressure from 50 kPa to 100 kPa, so that further 
precision on the neutron lifetime would be possible. 



Kossakowski et al. 

This work 

■Benv 

80 cps 

7.7 cps 

Bbeam 

20 cps 

2.0 cps 

£carbon 

11% 

0.1% 

Unc. on statistics 

3.1% 

1.0% 

Unc. on separation between iV3 He and Np 

0.6% 

0.6% 

Unc. on subtraction of B beam 

0.9% 

< 0.9% 

Unc. on correction of ep 

< 0.1% 

< 0.1% 

Unc. on correction of e backscatter 

0.5% 

0.5% 

Unc. on correction of e carbon 

1.0% 

< 0.1% 

Unc. on p 3 He 

0.5% 

0.3% 

Unc. on cr3 He n n 

0.1% 

0.1% 

Total uncertainty 

3.5% 

< 1.6% 


Table 4: Comparison of B env , Bbeam , e car bon and neutron lifetime uncertainties between 
the experiment conducted by Kossakowski et al. and the present experiment, assuming 
300-kW beam power at J-PARC and data acquisition time of 34 hours. 
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Appendix A. Modeling the saturation of multiplication 

Suppose that we have a single-wire drift tube with sense-wire radius, a, 
and tube radius, b. We take a cylindrical coordinate system with the z-axis 
aligned in the wire direction and consider a multiplication process problem 
in which a dense bunch of ionized electrons drift to a narrow spot on the 
wire. Assume that the process is sufficiently fast, i.e., the ions created by 
avalanche multiplication do not move away from their created position close 
to the wire in the time being considered. Also, we approximate the model by 
assuming that every ionized electron multiplication process is instantaneous 
and sequential in time; an ionized electron drifting to the wire only looks 
at the electrostatic held around the wire, which was formed by all previous 
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multiplications. The space charge distribution around the sense wire at the 
time that n ionized electrons have already multiplied is labeled N n {r). From 
Gauss’s law, the electric field, E n (r), is written as 

pr pd 

2irreE n (r) = ecrd / dr' / dz'2nr'N n (r'), (A.l) 

J a JO 

where e, e, cr, and d denote the dielectric constant, elementary charge, line 
charge density on the wire, and system dimension in the ^-direction, respec¬ 
tively. The boundary condition is 

drE n (r) = V 0 , (A.2) 

where Vo is the high voltage applied to the sense wire. The form of N n (r ) 
is not known precisely but, if we assume that the shape of the distribution 
does not change drastically during the multiplication process, we can write 

N n(r) = '^( r )’ ( A - 3 ) 

where Gj is the multiplication gain for the j-th ionized electron and p(r) is 
the average space-charge distribution, normalized by 

dr p(r) = 1 . (A.4) 




Then, the electric field is formally expressed as 

' n —1 


EJr) = 


v q , _i_ / vr 

rln(b/a) ead J 

/ P(r') 


x 


- / dr' + p( r ) 

l ) J a r 


(A.5) 


[ In (b/a) 

where the first term is the electric field with no space charges, the second 
term is the “screening effect” of the electric field due to the space charges, 
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and the third term is the Coulomb held created by the space charges. We 
do not calculate the multiplication gain, G n , for the given E n (r ) explicitly. 
Alternately, we suppose a real value, E n , which characterizes the scale of 
E n (r). Under the assumption that p{r) does not change, we can define E n 
as the value E n (r = a). For example 


E n = E n (r = a) = 


Vn 


eC 


' n— 1 


c = 


In(6/a) J a 


aln(6/a) ead Q 
' d(r') 

dr' — - - p(a) = const 


E G i b 


Then, we employ an empirical relation between G n and E n , 
widely applicable to wire chambers 


(A.6) 

(A.7) 
which is 


G n = exp (A + BE,,) , 


(A.8) 


where A and B are arbitrary constants used to characterize the multiplica¬ 
tion. From Eq. (IA.6[) and (1A.8D . if we take the difference, E n — E n _ i, we can 
deduce a recurrence relation for the gain such that 

G n = G n -1 exp ^—gad ' (A.9) 

The above recurrence relation is approximately solved in the case of aG o = 
eBC/(ead ) • G 0 1 to 


G r , 


(A.10) 


Go 1 T oiGqti 

By summing over all N ionized electrons and normalizing by the total multi¬ 
plication, assuming no space-charge effects, the degree of degradation of the 
multiplication caused by the space-charge effect is expressed as 


s(G 0 ) = 



ln(l + aNG 0 ) 
aNGo 
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ln(l + /A^Go) 

fAEG 0 


(A.ll) 


where A E is the amount of energy deposit observed by the wire, which is 
proportional to N. The coefficient, /, is hence defined as f = aN/ AE. We 
call / the saturation parameter. Eq. flA.lip indicates that / is a measurable 
value, provided A E has been determined by scanning Go and relatively com¬ 
paring the variance of s(Go). We performed the measurement of s(Go) for 
our TPC with configuring a 241 Am source collimated to fixed incident angles. 
Figure IA.10I shows an example of the scan of s(Go) by changing the sense 
wire’s HV. The gain Go at a given sense wire’s HV was monitored using 
a 55 Fe 5.9 keV X-ray source. The common absolute normalization for the 
set of the HV scan and f AE are fitted simultaneously following Eq. flA.lip 
(the vertical scale of Figure IA.10I is already normalized using the fitted ab¬ 
solute scale so that s(Go) becomes 1 at Go —>• 0). The deposit energy on the 
wire AE is undetermined from the measurement, and it was estimated us¬ 
ing Geant4 simulation by comparing the energy deposit on different MWPC 
wires to extract the saturation coefficient /. Once / is obtained, it is pos¬ 
sible to calculate the space-charge effect as a function of the primary gain, 
G 0 . Eq. flA.lip does not explicitly depend on the precise structure of the 
single-wire drift tube (e.g., parameters such as the wire radius), and we can 
obtain the same equation for MWPCs; the precise structure related to the 
saturation is all included in the saturation parameter, /. 

Let us then consider determining / in the case of our TPC, for a dense 
track. We use the same coordinate system defined in Section [2] It is obvious 
that / depends on the effective length of the wire, which is affected by the 
space-charge effect. Provided the drifting of positive ions away from the sense 
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Figure A.10: Measurement of the TPC’s saturation s(Go) as a function of gain for a 
certain energy deposit on a sense wire given by the 241 Am source injection. The vertical 
scale was normalized to s(Go) —> 1 at the limit of Go —» 0 by fitting the data points with 
Eq. [ATl] multiplied by an arbitrary scale. The product of the saturation parameter / and 
the energy deposit A E is obtained by the fitting. 


wire is negligible, we can neglect the time dependence, and the system is then 
two-dimensional. Thus, / should be a function of the angle of the particle 
track, 6, projected onto the xz -plane with respect to the wire. Qualitatively, 
we can say that / takes the largest value at 9 = 7r/2 and the smallest at 
9 = 0. Also, / should be even for 9 < 7t/2 and 9 > tt / 2 , and it is smooth 
and at 9 = tt/2. The simplest form of f{9) will be written as a function of 


m 


fn 


y/l + (/CCOt 9) 2 ’ 


(A-12) 


where / max is the maximum saturation parameter at 9 = tt/2 and k is the 
strength of the d-dependence. These parameters can be obtained from the 
actual measurement by using a 241 Am source, as shown in Figure lA.lll Math¬ 
ematically, the energy deposit density on the wire diverges at 9 = tt/2 but, in 
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reality, it remains finite due to multiplication dispersion in the z-direct ion, as 
well as the diffusion of ionized electrons before multiplication. At the 6 —» 0 
limit, / —>• 0 and s(Go) —> 1 and no space-charge effect exists. 



0 (deg) 


Figure A.11: Saturation parameter, /, as a function of the particle track angle, 0, projected 
onto the rrz-plane with respect to the wire. The parameters, / max and k, are obtained 
from the measurements indicated by circles. 
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